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Abstract. Noncollinear magnetic interfaces introduced in superconductor 
(SC)/ferromagnet/SC heterostructures allow for spin-flipping processes and 
are able to generate equal-spin spin-triplet pairing correlations within the 
ferromagnetic region. This leads to the occurrence of the so-called long-range 
proximity effect. Particular examples of noncollinear magnetic interfaces include 
Bloch and Neel domain walls. Here, we present results for heterostructures 
containing Bloch and Neel domain walls based on self-consistent solutions of the 
spin-dependent Bogoliubov—de Gennes equations in the clean limit. In particular, 
we investigate the thickness dependence of Bloch and Neel domain walls on 
induced spin-triplet pairing correlations and compare with other experimental 
and theoretical results, including conical magnetic layers as noncollinear magnetic 
interfaces. It is shown that both, Bloch and Neel domain walls lead to the 
generation of unequal-spin spin-triplet pairing correlations of similar strength as 
for conical magnetic layers. However, for the particular heterostructure geometries 
investigated, only Bloch domain walls lead to the generation of equal-spin spin- 
triplet pairing correlations. They are stronger than those generated by an 
equivalent thickness of conical magnetic layers. In order for Neel domain walls to 
induce equal-spin spin-triplet pairing correlations, they have to be oriented such 
that the noncollinearity appears within the plane parallel to the interface region. 
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1. Introduction 

Bringing a superconductor (SC) and a ferromagnetic (FM) material in close proximity 
to each other results in drastic changes for the spin-singlet Cooper pairs. While 
the Pauli principle requires the spins to orient parallel in the FM region, the spin- 
singlet Cooper pair electrons tend to align antiparallel in the SC region. This leads 
to the following consequences: i) Different Fermi velocities in the spin-up and spin- 
down channel of the electrons due to the exchange field lead to a centre of mass 
modulation and an oscillating behaviour of the superconducting correlations within 
the FM region Ells], namely the FFLO oscillations EE]. Unequal-spin spin-triplet 
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correlations {S^ = 0) are also generated. Both of these correlations are oscillating 
and suppressed in the FM region, and are essentially short-range, ii) Following a 
theoretical suggestion by Bergeret et al. [5] it should also be possible to induce equal- 
spin spin-triplet Cooper pairs which are compatible with and unaffected by the FM 
exchange field, thus leading to much larger penetration depths. This phenomenon 
requires spin-flip processes at the interface and is called the long-range proximity 
effect. 

Since this first theoretical prediction equal-spin spin-triplet pairing correlations 
have attracted a lot of attention and are reviewed by Buzdin [7], Bergeret et al. [5], 
and Tanaka et al. [3. Several multilayer setups have been suggested experimentally 
and theoretically to shed some light on the specific interface properties necessary to 
generate equal-spin spin-triplet pairing correlations within such structures. These 
setups typically involve noncollinear magnetic structures within the interface region 
as provided by different FM layers [TUI [H] ; helical (or conical) magnetic layers [TUI 
II1[I1[IS1[IU], or through specific interface potentials [I7||I1II11[1U]. Moreover, also 
the effects of Bloch [Ml [151 HI] and Neel [Ml [151 HU HI] domain walls, as well as 
interfacial spin-orbit coupling [231124] as source for generating equal-spin spin-triplet 
pairing correlations have been investigated. 

On the theoretical side, equal-spin spin-triplet pairing correlations have been 
investigated using Green’s function techniques based on solutions of the Usadel 
equations [HlTllTSlllT], and self-consistent solutions of the Bogoliubov—de Gennes 
(BdG) equations [161 Hi] • While the Green’s function techniques allow for an easier 
inclusion of scattering effects describing a shorter electron mean free path in the 
diffusive limit, these effects can in principle be incorporated into the BdG equations 
as well via additional impurity scattering potentials. However, in the present work we 
are solely concerned with the special case of odd-frequency triplet pairing correlations, 
which are not affected by nonmagnetic impurities [8]. The inclusion of additional 
impurity scattering potentials would require further investigations. 

The main aim of the paper is to assess the effectiveness of Bloch and Neel domain 
walls in comparison to conical magnetic layers to generate long-range proximity 
effects. This is motivated by the potential use of magnetic domain walls incorporated 
in experimental heterostructures to allow for a switching of the generated long- 
range spin-triplet Josephson current. Based on an established heterostructure setup 
successfully used in previous experimental and theoretical works we introduce Bloch 
and Neel domain walls in the interface regions and investigate its suitability to generate 
equal-spin spin-triplet pairing correlations. These results will be compared to the 
corresponding results with interface regions containing conical magnetic layers of 
varying thicknesses. 

This paper is organised as follow. Section [2] gives an overview over the theoretical 
methods. This includes a description of the microscopic BdG equations in section [2Tl 
and the definition of the spin-triplet pairing correlations in section 12.21 The results 
for heterostructures containing Bloch and Neel domain walls of varying thicknesses 
’T'Bioch and UNeei are presented in section 13.11 and section 13.21 respectively. Section |4| 
provides a summary and an outlook. 
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2. Theoretical background 

2.1. BdG equations and heterostructure setup 


The results presented in this paper are obtained by means of the microscopic BdG 
equations, which have been solved self-consistently in the clean limit. In the most 
general spin-dependent case and incorporating an expression for arbitrary exchange 
fields h needed later to describe the influence of Bloch and Neel domain walls, the 
BdG equations read [ISlIlilllS] 
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where e„, Una, and Vna denote eigenvalues, and quasiparticle and quasihole amplitudes 
for spin tr, respectively. Following simplifications introduced earlier [M 123 [2H] the 
tight-binding Hamiltonian TLq in © becomes effectively one-dimensional and reads 

TLq = —t 'y ^ ^cj^Cn+l + ^ ^ {Sn ~ M) . (2) 


Therein, ct and c„ denote electronic creation and destruction operators at multilayer 
index n, respectively. In order to compare the present results with those obtained 
previously [161 IMl |30] we choose the next-nearest neighbour hopping parameter t = 1 
and set the energy scales via the chemical potential (Fermi energy) p, = 0. 

Balian and Werthamer ISIISI] introduced a way to rewrite the general form of 
the pairing matrix appearing in © as 


(3) 


A-j-i- \ /A,^ f + idy A + dz 

Here, denotes a 2 x 2 matrix. Making additional use of the Pauli matrices a, 
the superconducting order parameter is now written as a singlet (scalar) part A and a 
triplet (vector) part d, respectively. This work is restricted to s-wave superconductors, 
and the pairing potential entering ([1]) simplifies to a scalar quantity A, which fulfills 
the self-consistency condition 

A(r) = ^ (^^"t(i’Xi(i’)[l - /(en)] + Uni{r)v*^{r)f{en)) ■ (4) 


The sum is evaluated only over positive eigenvalues e„ and /(Cn) denotes the Fermi 
distribution function evaluated as a step function for zero temperature. 

The heterostructure is set up as generally shown in figure [11(a) with a lattice 
constant a = 1 and a strength of the exchange field ho = 0.1. The effective 
superconducting coupling parameter g{r) equals 1 in the nsc = 250 left and right 
layers of spin-singlet s-wave superconductor and vanishes elsewhere. Between the 
s-wave superconductors and npM = 100 layers of the ferromagnetic middle layer 
we introduce conical magnetic layers (figure [IKa)), and Bloch (figure HKc)) and Neel 
(figure [Tl(d)) domain walls of varying thicknesses tt-cm, ?^Bioch, and riNeei, respectively. 
The conical magnetic layers are arranged as in our previous works [ISllao] and are 
chosen to represent the conical magnet Holmium. The maximum number of conical 
magnetic layers considered here correspond to one full turn of the magnetisation 
around the cone. The Holmium opening angle a = 80° and turning (or pitch) angle 
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(a) 







Figure 1. Structural setup of the heterostructures investigated in this work: (a) 
General setup including nsc layers of s-wave superconductor to the left, ncM 
layers of a conical magnet, npM layers of ferromagnetic middle layer, and ncM 
conical magnet and nsc s-wave superconductor layers to the right, (b) Employed 
cordinate system and definition of the opening angle a and the turning (or pitch) 
angle ^ of the conical magnetic structure, (c) Orientation of the Bloch domain 
wall of nBioch layers between the left s-wave superconductor and the middle 
ferromagnetic region, (d) Orientation of the Neel domain wall of n^^el layers 
between the left s-wave superconductor and the middle ferromagnetic region. 
Note that (c) and (d) only show the left interface region for clarity. 


P = 30° are measured from the positive y axis towards the positive ^ axis and from 
the positive ^ axis towards the positive x axis, respectively (figure [T{b)). The Bloch 
and Neel domain walls are arranged such that there is a ferromagnetic coupling to the 
middle ferromagnetic region of the heterostructure and that a C 2 symmetry along the 
central z axis is retained. 

2.2. (Triplet) pairing correlations 

The specific superconducting pairing correlation between spins a and /3 we are 
interested in here is evaluated as on-site interaction for times t = t and t' = 0 

/c/3(r,r,0) = r)4';3(r,0)) . (5) 

Therein, r) denotes the many-body field operator for spin a at time r, and time- 

dependence is introduced through the Heisenberg equation of motion. The particular 
pairing correlation of (O is local in space and leads to vanishing triplet contributions 
for r = 0, in accordance with the Pauli principle |33| . However, finite values of t give 
rise to nonvanishing pairing correlations, being an example of odd-frequency triplet 
pairing [5]. Substituting the field operators valid for our setup and phase convention 
the spin-dependent triplet pairing correlations read Hi [30] 

fo{y,T) = (u„-f(?/X^(2/) + u„;(2/X-,-(y))Cn(r) 

1 1 ” ’ 

fi(.y,r) = -^{fniy^T) - fuiy^T)) = 2N - wn;(2/)X(2/))Cn(r) 
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depending on position y within the heterostructure and time parameter r (set to 
T = 10 in the present work), and with Cn('r) given by 

C«(t) = cos(e„T) - zsin(e„T)(l - 2/(e„)) . (7) 

Based on the rewritten form of the pairing matrix in ([3]) and defining A as the triplet 
pairing matrix for an ordinary spin-triplet superconductor, the product AA^ can be 
written as 

AA'^ = |dp(7o + f (d X d*) a. (8) 

Therein, the magnitude of the d-vector denotes the gap function independently from 
the underlying coordinate system, and id x d* is a measure of the spin magnetic 
moment, respectively. In the present work there is no pairing interaction considered 
in the triplet channel, i.e., g{r) = 1 in (jdj). However, triplet pairing correlations 
are induced which are given by the so-called triplet pair correlation function matrix 
instead, which can be written similarly to the pairing matrix A as |341130] 

With © the analogue to (|5]) can be written as 

fP = \i\^ao+i{fxr)a. (10) 

Analogously to |d| and d x d* we find |f| and f x f*. Again, they can be rewritten in 
terms of the f-vector components which depend on the spin-dependent triplet pairing 
correlations as 

/x = 2 (“/n + fu) 

fy = -lifn + fu) ■ ( 11 ) 

/z = 2 ifn + At) 

3. Results and Discussion 

3.1. Heterostructure including Bloch domain walls 

The results presented in this section are obtained for heterostructures containing Bloch 
domain walls of increasing thickness riBioch to either side of the ferromagnetic middle 
layer. The Bloch domain walls are oriented such that the magnetic moments at the 
interface align parallel to the ferromagnetic middle layers (figure [IJc)) and that a C 2 
symmetry along the central z axis is retained. Recently, heterostructures containing 
conical magnetic layers of different opening and turning angles leading to different 
symmetries have been investigated [HI [30] and could be explained with the properties 
of the f-vector introduced in (|6|). The fundamental pairing correlations evaluated 
according to (jS]) are shown in figure [5] and figure |3I for heterostructures containing 
conical magnetic layers ncM, Bloch domain walls riBioch, and Neel domain walls riNeei, 
respectively. In particular, the middle panels of figure [2] show the magnitude of the 
unequal-spin spin-triplet pairing correlations fo for increasing thickness of the Bloch 
domain wall riBioch in comparison to conical magnetic layers ncM (upper panels) and 
Neel domain wall riNeei (lower panels), respectively. The full data sets of the left panels 
are shown as top views in the respective right panels. 
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Figure 2. Magnitude of the unequal-spin spin-triplet pairing correlations fo for 
heterostructures containing conical magnetic layers riQyi (upper panels), Bloch 
domain walls riBioch (middle panels), and Neel domain walls (lower panels), 

respectively. The left and right panels both show the full data sets and the right 
panels show a top view of the data. 


Looking at the oscillations within the ferromagnetic middle region of the 
heterostructure, one notices strong oscillations of the unequal-spin spin-triplet pairing 
correlations fo- This is very similar to heterostructures containing conical magnetic 
layers (shown in the upper panels). These are essentially the FFLO oscillations 
discussed above. The magnitude is almost equally strong for thin Bloch domain walls, 
but decays to smaller values for increasing Bloch domain wall thickness compared to 
the conical magnetic layers. Flowever, within the interface region the magnitude of 
fo is increasing linearly with Bloch domain wall thickness and clearly surpasses the 
strength in case of the conical magnetic layers. Looking now at figure [3] showing the 
respective results for the magnitude of the equal-spin spin-triplet pairing correlations 
fi for heterostructures containing conical magnetic layers ncM (upper panels) and 
Bloch domain walls riBioch (lower panels) one immediately notices larger differences. 
The magnitude of fi rises linearly and is similarly strong for smaller values of Bloch 
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Figure 3. Magnitude of the equal-spin spin-triplet pairing correlations fi 
for heterostructures containing conical magnetic layers ncM (upper panels), 
and Bloch domain walls naioch (lower panels), respectively. Note that the fi 
contribution from heterostructures containing Neel domain walls is zero for the 
chosen geometry (as discussed in the text). The left and right panels both show 
the full data sets and the right panels show a top view of the data. 


domain wall thickness neioch compared to the same thickness of conical magnetic 
layers ncM- In contrary to the conical magnetic layers, for which fi starts to 
oscillate and decay with a further increase of layer thickness again, the equal-spin 
spin-triplet pairing correlations fi for thicker Bloch domain walls increase further and 
reach saturation for the maximum thickness considered here. Moreover, the overall 
magnitude of fi for heterostructures containing Bloch domain walls clearly exceeds 
those for heterostructures containing conical magnetic layers. Finally, a measure for 
the induced spin magnetic moment provided by the quantity f x f* of cni) and m is 
shown in figure |4j where again the contributions from heterostructures containing 
conical magnetic layers ncM, and Bloch riBioch and Neel riNeei domain walls are 
shown in the upper, middle, and lower panels, respectively. One immediately notices 
that spin magnetisation only occurs in the respective interface regions and that no 
magnetisation leaks into the SC regions of the heterostructures. The overall rise in 
spin magnetisation is stronger for the Bloch domain walls and saturates earlier than 
in the conical magnetic layers. The latter also show stronger oscillations depending 
on the thickness uqm compared to the Bloch domain wall thickness riBioch- 

Summarising this we find a similar influence on fp for heterostructures containing 
conical magnetic layers and Bloch domain walls within the ferromagnetic middle layer, 
and a stronger influence within the Bloch domain walls for increasing thickness nsioch- 
The magnitude of fi for heterostructures containing Bloch domain walls riBioch clearly 
exceeds those for conical magnetic layers and saturates for the maximum thickness 
of UBioch considered here. The measure of induced spin magnetisation provided by 
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Figure 4. Magnitude of f X f* for heterostructures containing conical magnetic 
layers ncM (upper panels), Bloch domain walls riBioch (middle panels), and Neel 
domain walls (lower panels), respectively. The left and right panels both 

show the full data sets and the right panels show a top view of the data. 


f X f* rises faster and saturates earlier in heterostructures containing Bloch domain 
walls compared to conical magnetic layers. 

3.2. Heterostructure including Neel domain walls 

Now we focus on results obtained for heterostructures containing Neel domain walls 
of increasing thickness riNeei to either side of the ferromagnetic middle layer. They 
are oriented such that magnetic moments of the Neel domain wall align parallel to the 
ferromagnetic middle layers at the interface (figure [TJd)) and are oriented in the yz 
plane. 

Starting the discussion again with the unequal-spin spin-triplet pairing 
correlations fp shown for Neel domain walls of varying thickness riNeei in the lower 
panels of figure one notices that they are nearly indistinguishable to the respective 
contributions from the Bloch domain walls shown in the middle panels of figure [2l 
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It seems that the orientation of the noncollinear magnetic moments with respect to 
the interface plane does not influence the induced unequal-spin spin-triplet pairing 
correlations. 

This changes drastically when looking at the equal-spin spin-triplet pairing 
correlations where we find no contribution stemming from heterostructures containing 
Neel domain walls (within the particular setup investigated in the present work). This 
phenomenon has been discussed recently by Alidoust et al. [14] based on the properties 
of the triplet anomalous Green’s function f used in their approach. It turns out that for 
Neel domain walls to successively induce equal-spin spin-triplet pairing correlations fi 
the respective noncollinear magnetic structure has to occur in the plane parallel to the 
interface. In the present work the noncollinearity of the Neel domain wall is oriented 
perpendicular to the interface plane and hence does not contribute towards equal-spin 
spin-triplet pairing correlations. 

Focusing now at the properties of f x f* shown in figure S] one notices again 
that the contributions from heterostructures containing Bloch nsioch (middle panels) 
and Neel riNeei domain walls (lower panels) are again indistinguishable. The spin 
magnetisation is only present in the interface regions and does not leak into the SC 
regions of the heterostructures. Again, the overall rise in spin magnetisation is stronger 
for the Neel domain walls and saturates earlier than in the conical magnetic layers 
(upper panels). The latter also show stronger oscillations depending on the thickness 
ncM compared to the Neel domain wall thickness riNeei- 

In summary, for the heterostructures containing Neel domain walls we find 
indistinguishable results for the unequal-spin spin-triplet pairing correlations fo and 
the measure of the spin magnetisation f x f* compared to heterostructures containing 
Bloch domain walls as discussed in section o However, for the setup chosen for the 
Neel domain walls we find no induced equal-spin spin-triplet pairing correlations fi at 
all. 

4. Summary and Outlook 

Here we presented a detailed analysis of spin-triplet pairing correlations occurring 
in SC/FM/SC heterostructures incorporating different types of noncollinear interface 
regions of varying thicknesses. Particular interfaces included conical magnetic layers, 
and Bloch and Neel domain walls. It has been shown that for the unequal-spin spin- 
triplet pairing correlations fg both Bloch and Neel domain walls exceed the efficiency 
provided by conical magnetic layers at the interface. For the equal-spin spin-triplet 
pairing correlations fi only the Bloch domain walls show considerable improvement 
over the conical magnetic layer setup, whereas the specific orientation of the Neel 
domain wall (noncollinear magnetic moments perpendicular to the interface plane) 
prevented the occurrence of equal-spin spin-triplet pairing correlations in that case. 
Different orientations of Neel domain walls with respect to the interface plane deserve 
further investigations. Finally, the dependence of the spin magnetisation f x f* on the 
thickness of both, Bloch and Neel domain walls, appears to be smoother compared to 
the conical magnetic layer setup, but in all cases this stays restricted to the interface 
region alone. 
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